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Abstract 

We introduce a new Monte Carlo method for pure gauge theories. It is 
not intended for use with dynamical fermions. It belongs to the class of Local 
Hybrid Monte Carlo (LHMC) algorithms, which make use of the locality of the 
action by updating individual sites or links by following a classical mechanics 
trajectory in fictitious time. We choose to update a one-parameter subgroup 
of the gauge field on each link of the lattice, and the classical trajectory can 
be found in closed form in terms of elliptic functions for this case. We show 
that this gives an overrelaxation algorithm with a tunable parameter which, 
unlike some previous methods, does not require the numerical integration of 
the equations of motion. 



*To appear in the Proceedings of Lattice '93, Dallas, TX; October 1993. 
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1 LHMC FOR FREE FIELD THEORY 

Let us first quickly review the Local Hybrid Monte Carlo (LHMC) and Adler's over- 
relaxation (AOR) algorithms [1] for free field theory. See [2] for further references. 
Consider the Gaussian model defined by the free field action 



A single-site AOR update with parameter ( replaces by 



(j) [x) = (1 - Cmx) + ^ + V, 



where uj^ = 2D + m^, J- = J2\x-y\=i 4'{.y)i ?/ is a Gaussian random number. 
The lattice may be updated using a checkerboard scheme. 

• For C = 1 this is the heatbath (HB) algorithm, because (j)'{x) does not depend 
upon (j){x). 

The exponential autocorrelation time is 

D 

vn? + 

which corresponds to z = 2. 



= -T—r, {f - 0) 



• If we adjust ^ so as to minimize the autocorrelation time we find that 



C = 2 - ^ + 0{m^), 



which gives z = 1. 

Consider the LHMC algorithm applied to the Gaussian model: We introduce the 
Hamiltonian H = |p^ + S'(0) on "fictitious" phase space. The corresponding equation 
of motion are 

4>x = —i^^4>x + ^, 

whose solution in terms of the Gaussian distributed random initial momentum pg. and 
the initial field value (j)g. is 

1 — cos cut Pg. 
(pJt) = (pg. cos u)t H T H sm u)t, 

which is exactly the AOR update considered before if we identify ^ = 1 — cos Lot and 
Vx = V- 

If we use this exact solution of the equations of motion to generate candidate 
configurations for LHMC then the acceptance rate will be unity. 
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2 LHMC FOR PURE GAUGE THEORY 



We now turn to the case of pure gauge theory. The Wilson action is 



S{Ue) = -/3Re 



TtV 



where the matrix V = J2ved*i U-p. 

Our update will be Ui gUi where 5^ is a group element. 

There is no reason why all the degrees of freedom on a link need to be updated 
simultaneously; we are free to update just a single degree of freedom at a time as long 
as we can guarantee ergodicity. Therefore we shall restrict g to lie in a one- dimensional 
subgroup. 

For U in the fundamental representation of SU{N) or SO{N) any group element 
can be written as a conjugate of 



-ie 



or 



cos 6 sin 6 
— sin 6 cos 6 



respectively, so updating with just this subgroup will be ergodic if random gauge 
transformations are carried out from time to time. In other words, selecting an 
element from these subgroups is equivalent to choosing a completely general group 
element up to gauge transformations. 
If we introduce the parameters 



ae 



{ Vii + V*^ for SU{N) 

i{Vr2-V2i) ioT SO{N) 



then we may write the action in the simple form 

s{e) = ^{1 - cos(^ + s)}. 



As usual, we introduce a fictitious momentum p and a Hamiltonian 

H{p,e) 



where Vq = 2/3a/N is the maximum potential energy of the "pendulum". 

Our LHMC algorithm consists of repeatedly applying the following steps in any 
order: 



• Refresh the momentum p from a Gaussian heatbath; 
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• Integrate the equations of motion for a time t[E) and then reverse the sign of 
the final momentum; 

• Perform a random gauge transformation. 

The only subtlety is that we can allow the trajectory length r depend upon the 
constants of motion (specifically the energy E) because for an exact solution of the 
equations of motion this still satisfies detailed balance. 

We now have to solve the equations of motion, which are 

e = -VoS^n ^-^j cos ^-^j 

subject to the initial conditions ^(0) = and p(0) = po- 

It is convenient to introduce a new variable y = sin and a rescaled time 

dy 



^ = ^(l-y2)(l-fcV) 



u = tJE l2, whence 



where the parameter k = ^Vo/ E. The solution of this equation is the Jacobian elliptic 
function 

y = sn('iisgnpo + Uo,k) 

where snuo = sin |. 

Remember that we can choose any value u[k) for the trajectory length. One 
particularly interesting choice is to take the trajectory length to be a fixed fraction 
of half the period of an orbit, u, which corresponds to overrelaxation; 

u{k) = ^1 - 

where ^ is (proportional to) the longest relevant correlation length of the gauge theory. 

A pendulum exhibits two different types of motion depending on the value of the 
parameter k (for k = 1 the solution is a Sine-Gordon instanton) 



2K{k) iik <1, E >Vo 
iK{^ iik>l,E<Vo, 

where K is the complete elliptic integral 

dy 



•^0 ^/ 1 



/(l-y2)(l-fcV) 

With this choice of trajectory length the value of 6 selected may be written as 

- / usgnpo \ 
a = Z am I j j , 

where the amplitude a.m.u = sin~^ snu. 
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H(q,p) 




Figure 1: The Hamiltonian as a function on phase space (with q = 6) with the 
corresponding classical trajectories. 
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Figure 2: Preliminary results for two-dimensional SU{2) on an 128^ lattice at /3 = 
20.0. Convergence for various values of m = 1/^ and for hot and cold starting 
configurations is demonstrated. Studies of the dynamical critical behaviour are in 
progress. 
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3 CONCLUSIONS 



We may hope that the LHMC algorithm, with optimally chosen trajectory lengths, 
will exhibit a dynamical critical exponent z ^ 1, ii the results from free field theory 
prove to be applicable to pure gauge theories. In this regard it is similar to the Hybrid 
Overrelaxation algorithm analysed by Wolff. 

The advantages of the LHMC method as compared to previous overrelaxation 
methods for pure gauge theories is that it allows a true tunable overrelaxation pa- 
rameter without the need for explicit numerical integration. 
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